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Abstract 

We report measurements of the exclusive electroproduction of K~^K and K^TP final states from a 
proton target using the CLAS detector at the Thomas Jefferson National Accelerator Facility. The 
separated structure functions ctt, ct^, gtt^ and glt were extracted from the and e-dependent 
differential cross sections taken with electron beam energies of 2.567, 4.056, and 4.247 GeV. This 
analysis represents the first (Tl/(Tt separation with the CLAS detector, and the first measurement 
of the kaon electroproduction structure functions away from parallel kinematics. The data span 
a broad range of momentum transfers from 0.5 < < 2.8 GeV^ and invariant energy from 
1.6 <W < 2.4 GeV, while spanning nearly the full center-of-mass angular range of the kaon. The 
separated structure functions reveal clear differences between the production dynamics for the A 
and hyperons. These results provide an unprecedented data sample with which to constrain 
current and future models for the associated production of strangeness, which will allow for a better 
understanding of the underlying resonant and non-resonant contributions to hyperon production. 

PACS numbers: 13.40.-f, 13.60.Rj, 13.85.Fb, 14.20.Jn, 14.40.Aq 
Keywords: CLAS, kaon electroproduction, structure functions, hyperons 
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I. INTRODUCTION 



A necessary step toward understanding the structure and dynamics of strongly interacting 
matter is to fully understand the spectrum of excited states of the nucleon. This excitation 
spectrum is a direct reflection of its underlying substructure. Understanding nucleon reso- 
nance excitation, and hadro-production in general, continues to provide a serious challenge 
to hadronic physics due to the non-perturbative nature of the theory of strong interactions. 
Quantum Chromodynamics (QCD), at these energies. Because of this, a number of approx- 
imations to QCD have been developed to understand baryon resonance decays. One such 
approach is a class of semi-relativized symmetric quark models that invoke massive 

constituent quarks. These models typically predict many more nucleonic states than have 
been found experimentally. A possible explanation to this so-called "missing resonance" 
problem is that these nucleon resonances may have a relatively weak coupling to the pion- 
nucleon states through which many searches have been performed, and may, in fact, couple 
to other final states such as multi-pion or strangeness channels. In this work we provide 
an extensive set of data that may be used to search for these hidden states in strangeness 
electroproduction reactions. These data then provide for a complementary way in which to 
view the baryon resonance spectrum, as some of the "missing" states might be only "hidden" 
when studied in particular reactions. It could also be the case that some dynamical aspect 
of hadronic structure is acting to restrict the quark model spectrum of states to the more 

n 

limited set established by existing data p|. 

Beyond different coupling constants relative to single-pion production (e.g. Qkyn vs. 
Q-kMn)-, the study of the exclusive production of K^K and K^TP final states has other 
advantages in the search for missing resonances. The higher masses of the kaon and hyperons, 
compared to pionic final states, kinematically favor a two-body decay mode for resonances 
with masses near 2 GeV, a situation that is experimentally advantageous. New information 
is also provided by comparing K^K to K^TP . Note that although the two ground-state 
hyperons have the same valence quark structure (uds), they differ in isospin, such that 
intermediate N* resonances can decay strongly to K^A final states, but intermediate A* 
states cannot. Because K~^I1^ final states can have contributions from both A^* and A* 
states, the hyperon final state selection constitutes an isospin filter. 

Electro-excitation of the nucleon has served to quantify the structure of many excited 
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states that decay to single pions. These studies have been used to test models of the internal 
structure of the excitations using the interference structure functions and the behavior with 
for various regions of W . This experiment is the first to allow analogous investigations 
to start when the intermediate electro-excited nucleon resonances couple to kaon-hyperon 
final states. 

The search for missing resonances requires more than identifying features in the relevant 
mass spectrum. It also requires an iterative approach in which experimental measurements 
constrain the dynamics of various hadrodynamic models. The tuned models can in turn be 
used to interpret s-, t- and w-channel spectra in terms of the underlying resonances. As 
emphasized by Lee and Sato [4], QCD cannot be directly tested with A^* spectra without 
a model for the production dynamics. The key to constraining models and unraveling the 
contributing resonant and non-resonant diagrams that contribute to the dynamics, is to 
measure as many observables over as wide a kinematic range as possible. 

In this paper, we present measurements of the separated structure functions cif/, arr^ and 
a^T for exclusive electroproduction of K^K and K^TP final states for a range of momentum 
transfer from 0.5 to 2.8 GeV^ and invariant energy W from 1.6 to 2.4 GeV, while spanning 
the full center-of-mass angular range of the kaon. Our center-of-mass angular coverage is 
unprecedented. These are the first data published on exclusive KY electroproduction that 
extend beyond very forward kaon angles in the center of mass. At one value of Q^, Q^= 
1.0 GeV^, we were also able to separate the unpolarized structure function, au, into its 
components ax and ai using a traditional Rosenbluth separation and also an alternative 
e — $ Rosenbluth technique (where e is the transverse polarization of the virtual photon 
and $ is the angle between the electron and hadron planes). In this alternative method, we 
obtain the four structure functions in a single fit. This extensive data set should provide 
substantial constraints on the various hadrodynamic models (discussed in Section [IV|) . Due 
to the very large number of analysis bins encompassed by this work, only aportion of our 
available data is included here. The full set of our data is available in Ref. |5|. 

After a brief review of the relevant formalism in Section [Tll an overview of previous ex- 
perimental work in this area in Section IIIIl and a brief review of the current theoretical 
approaches in Section IIVI we present our measurements made using the CLAS detector in 
Hall B at Jefferson Laboratory (JLab) in Sections |V] through IVIII In Section IVIIIt the 
results are examined phenomenologically and compared to predictions from several models 
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that have not been "tuned" to this data set. Finally, we present our conclusions regarding 
the potential impact of these data in Section IIX[ Our conclusions regarding the s-channel 
baryon spectrum are, unfortunately, rather limited. Real progress on identifying heretofore 
"missing" resonances will only result from a judicious fitting of the theoretical and phe- 
nomenological models to these data and the remainder of the world's data on these final 
states. 

II. FORMALISM 

In kaon electroproduction a beam of electrons with four-momentum pe = {Ee,Pe) is 
incident upon a fixed proton target of mass Mp, and the outgoing scattered electron with 
momentum p^' = {Ef,i,Pf,i) and kaon with momentum px = {Ex,Pk) are measured. The 
cross section for the exclusive /^"'"-hyperon state is then differential in the scattered electron 
momentum and kaon direction. Under the assumption of single-photon exchange, where the 
photon has four-momentum q = pe — Pe' = {^^q)-, this can be re-expressed as the product 
of an equivalent flux of virtual photons and the 7*p center-of-mass (cm.) virtual photo- 
absorption cross section as: 

dE^idVte'dVt]^ ~ dn*j^' ^ ' 

where the virtual photon flux factor F depends upon only the electron scattering process. 
After integrating over the azimuthal angle of the scattered electron, the absorption cross 
section can be expressed in terms of the variables Q^, W , 9^, and $, where = —Q^ is 
the squared four-momentum of the virtual photon, W = ^JM'^ + IM^v — is the total 
hadronic energy in the cm. frame, Q\ is the cm. kaon angle relative to the virtual photon 
direction, and $ is the angle between the leptonic and hadronic production planes (see 
Fig. [1]). After introducing the appropriate Jacobian, the form of the cross section can be 
written as: 

da ^ day 

dQ^dWdn*j^ ~ "Id*^' 

where 
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(2) 



w 



a 



Ml 



(3) 



is the flux of virtual photons, 



e = I i + 2i|Ltan^ 



is the transverse polarization of the virtual photon, and 6'e' is the electron scattering angle 
in the laboratory frame. 

hadron reaction plane 




FIG. 1: Kinematics for K~^Y (where Y is either a A or TP) electroproduction defining the angles 
9*^ and <I> with respect to the center-of-mass reference frame. 



After summing over the polarizations of the initial and final state electrons and hadrons, 
the virtual photon cross section can be written as: 



Ot + to L + tOxT COS 2$ 



VjaiT cos $. 



(5) 



K 



In this expression, the cross section is decomposed into four structure functions, ur, ul, arr, 
and (JiTi which are, in general, functions of Q^, W , and ^^only. Note that this convention 
for the differential cross section is not used by all authors 6|]. 

Each of the structure functions is related to the coupling of the hadronic current to dif- 
ferent combinations of the transverse and longitudinal polarization of the virtual photon. 
(Jt = + <^A.) is the differential cross section contribution for unpolarized transverse 
virtual photons. In the limit 0, this term must approach the cross section for un- 

polarized real photons which only have transverse polarization, ol is the differential cross 
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section contribution for longitudinally polarized virtual photons, cttt and a^T represent in- 
terference contributions to the cross section. axT = — cr±) is due to the interference of 
transversely polarized virtual photons and is due to the interference of transversely and 
longitudinally polarized virtual photons. Here a\\ and a± are the cross sections for virtual 
photons having their electric vector parallel and perpendicular to the hadronic production 
plane, respectively. Note that the term axr in electroproduction is related to the linearly 
polarized photon beam asymmetry in photoproduction experiments, which is defined as 
Yj = — (Jxt/'^t- 

For the remainder of this paper we will refer to (Tjj = (TT + ecri as the "unseparated" part 
of the cross section. The further decomposition of these structure functions into response 
functions, which can then be expressed in terms of either complex amplitudes or as multipole 
expansions, is given, for example, in Ref. [tI. In this paper we will compare theory to the 
structure function terms introduced above. 



III. PREVIOUS EXPERIMENTAL WORK 

Hyperon electroproduction in the nucleon resonance region has remained largely unex- 
plored. Several low-statistics measurements were carried out in the 1970's at Cambridge, 
Cornell, and DESY, and focussed mainly on cross section measurements to explore differ- 
ences in the production dynamics between A and S*^ hyperons. The first experiment was 
performed at the Cambridge Electron Accelerator |^] using small- aperture spectrometers in 
kinematics spanning below 1.2 GeV^, W from 1.8 to 2.6 GeV, and forward kaon angles 
(^1^ < 28°). It was noted that the K~^A channel dominated the fC+S^ channel, with signs 
of a large longitudinal component in the K~^A channel. Subsequent results from Cornell 
for < 2.0 GeV^ and iy=2.15 and 2.67 GeV confirmed this observation. In this paper 
we show that K^A dominance only occurs at forward kaon angles, and the longitudinal 
strength is only important at forward angles and higher W. 



An experiment from DESY lOl] used a large aperture spark-chamber spectrometer to 
measure both reaction channels at higher W (1.9 < W < 2.8 GeV) and lower (0.1 < 
< 0.6 GeV^). That experiment managed the first and utt separations, albeit 
with large error bars, few data points, and considerable kinematic extrapolations to extract 
results at fixed values of W, Q^, and t = {q — Pk)'^- The results were consistent with zero for 
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these interference cross sections due to the large uncertainties. In this paper we show the 
first measurements of the structure functions with enough precision to determine non-zero 
interference terms. 

Other measurements made at Cornell were reported ll| for kaons produced at very small 
angles relative to the virtual photon {O^^ < 15°). The results for W from 2.15 to 3.1 GeV 
included improved measurements of the dependence of the differential cross sections over 
the range from 0.6 to 4.0 GeV^, showing that the K'^H^ cross section falls off much faster 
than the K~^A cross section. At that time this was explained by a vector meson dominance 
argument, or alternatively, as possible evidence for an isoscalar di-quark interaction that 
favors A production over production off the proton jl2|. Another survey experiment 
from DESY [13] at W=2.2 GeV and 0.06 < < 1.35 GeV^ that measured differential 
cross sections for the K~^A and K'^'E^ final states, confirmed the measured dependence 
of the cross sections, but with improved statistics. Our study of the dependence shows 
that this same behavior for S° production relative to A production also occurs for larger 
kaon angles. 

The first o"i/(TT separation via the Rosenbluth method was also made at Cornell 11| . 
suggesting that ctl for the K^A channel is large but not dominant at forward kaon angles, 
as previously surmised, while it is vanishing for the i^'+S" channel. At JLab two more 
recent results employing the Rosenbluth technique in parallel kinematics (0^=0°) have been 
completed to separate o"l and 0"^. The first result reported and cr^ for both the K~^A 
and final states using the small-aperture spectrometers in Hall C [l^. Results were 

extrapolated to W ~1.84 GeV for a range of from 0.52 to 2.00 GeV^. It showed, contrary 
to previous findings, that the ratio (Jl/(Jt for the TP is not very different in the forward 
direction than for the A over this range. The ratio for both hyperons is about 0.4, albeit 
with large uncertainties. The other existing measurement from JLab was performed in Hall 
A for W in the range from 1.8 to 2.14 GeV with Q"^ values of 1.9 and 2.35 GeV^ Q. This 
measurement was only for the K^A final state and showed that the ratio of (Tl/<^t was 
consistent with the Hall C result. 

In a previous CLAS publication using the same data presented in this work, the polariza- 
tion transfer from the virtual photon to the produced A hyperon was reported [IQ]. These ob- 
servables were expected theoretically to have strong sensitivity to the underlying resonance 
contributions. Surprisingly, they seemed to have only a modest dependence on W . The 
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CLAS polarization data were also analyzed to extract the ratio (Jl/<^t at 6*j^ = 0° [17||. The 
measured ratio was smaller than, but consistent with, that from the Hall C measurements, 
providing an important cross-check on the extraction of cr^ and cr^ from a measurement with 
very different systematics. 

In this paper we present data for ctl and ctt in similar kinematics that can be compared 
to these data. In addition we present the first available data for these structure functions 
for large 6^, away from parallel kinematics. Here, the longitudinal and transverse structure 
functions are extracted using the standard Rosenbluth technique, as well as by a simultane- 
ous €-$ fit to our different beam energy data sets. This analysis represents the first ctl/<^t 
separation using the CLAS spectrometer. 

In contrast to the sparse extant electroproduction data, there exist several high-quality 



photoproduction data sets. Recently, exclusive photoproduction o 



f K+A and ir+E° 



18 



inal 



19j and CLAS |20 



2l| 



states have been investigated with the large-acceptance SAPHIR 
detectors. High statistics total cross sections, differential cross sections, and induced polar- 
izations for the final state hyperons have been measured that span the full nucleon resonance 
region. In addition, high statistics measurements from CLAS of the beam-recoil hyperon 
polarization transfer have been completed for both the K^A and i^+S° final states 221], and 
)eam spin asymmetry measurements have been made at LEPS for both A and S hyperons 
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2J|. The W dependence of these data have been studied with the aim to understand 



the underlying s-channel A^* and A* contributions. Further information regarding interpre- 
tations of these data within different models is included in Section HVl 

Given this landscape of available data on the associated production of hyperons, it is clear 
that the majority of the existing electroproduction data, while spanning similar ranges of W 
and as our data, only provide information for very forward kaon scattering angles. This 
new data from CLAS represents a significant improvement in that it covers the full kaon 
scattering angular range, which will allow for an in-depth investigation of the contributing 
s-channel and u-channel diagrams in addition to the t-channel processes to these reactions. 
The new CLAS data also provides full azimuthal coverage, which enables the first signifi- 
cant data sample to study the interference structure functions. These structure functions 
provide new and unique information on interference between the underlying resonant and 
non-resonant amplitudes. In addition, CLAS has made significant contributions to the data 
base with the photoproduction cross sections and polarization observables that have been 
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published. This new set of electroproduction data allows the study of the production dynam- 
ics as a function of the mass of the virtual photon, which provides an exciting complement 
to the real photon data. 



IV. THEORETICAL MODELS 



At the medium energies used in this experiment, perturbative QCD is not yet capable of 
providing any analytical predictions for the differential cross sections or structure functions 
for kaon electroproduction. In order to understand the underlying physics, effective models 
must be employed that ultimately represent approximations to QCD. This paper compares 
the data to two different theoretical model approaches: hadrodynamic models and models 
based on Reggeon exchange. 



A. Hadrodynamic Models 

Hadrodynamic models provide a description of the reaction based upon hadronic degrees 
of freedom. In this approach, the strong interaction is modeled by an effective Lagrangian, 
which is constructed from tree-level Born and extended Born terms for intermediate states 
exchanged in the s, t, and u reaction channels as shown in Fig. O Each resonance has 
its own strong coupling constants and strong decay widths. A complete description of the 
physics processes requires taking into account all possible channels that could couple to the 
initial and final state measured, but the advantages of the tree-level approach are to limit 
complexity and to identify the dominant trends. In the one-channel, tree-level approach, 
several dozen parameters must be fixed by fitting to the data, since they are poorly known 
and not constrained from other sources. Identification of the important intermediate states 
or resonances is guided by existing data and quark model predictions. The coupling constants 
for each of the included resonances are extracted from global fits of the model calculations 
to the existing data base. It is common practice to use phenomenological form factors 



to account for the extension of the point-like interactions at the hadronic vertices |25 |. 
Different models typically have different prescriptions for restoring gauge invariance. The 
drawback of these models is the large number of exchanged hadrons that can contribute in 
the intermediate state of the reaction. Depending on which set of resonances is included, 
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very different conclusions about the strengths of the contributing diagrams may be reached. 
As stated in Section [H the models that are employed in this work have not been "tuned" 
to our data. It should also be stated that due to the nature of these models, they have a 
much higher interpretative power than predictive power. Therefore it will be the case that 
more definitive statements regarding the reaction dynamics and underlying resonant and 
background terms will only be possible after our data have been included in the model fits. 
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s-channel t-channel u-channel 

FIG. 2: Feynman diagrams representing s-channel nucleon {p, N* , A*) exchange (left), t-channel 
kaon {K,K*) exchange (middle), and u-channel hyperon {Y,Y*) exchange (right) that contribute 
to the reaction models. The vertex labels qmbb represent the strong coupling constants. 

Two different hadrodynamic models are employed in this work. The first is the model of 
Bennhold and Mart j26| (referred to here as BM) and the second is model B of Janssen et 



al. 



251] (referred to here as JB). In these models, the coupling strengths have been determined 



mainly by fits to existing 7p — * K^Y data (with some older electroproduction data included 
in some cases) by adding the non-resonant Born terms with a number of resonance terms 
in the s, t, and u reaction channels, leaving the coupling constants as free parameters. 
The coupling constants are required to respect the limits imposed by SU(3) allowing for a 
symmetry breaking at the level of about 20%. Both models have been compared against the 
existing photoproduction data from SAPHIR 19] and CLAS 20, 21] and provide a fair 
description of those results. The model parameters are not based on fits to any CLAS data. 
The specific resonances included within these calculations are listed in Table [H 

For K^K production, the BM model includes the Born terms as well as four baryon 
resonance contributions. Near threshold, the steep rise of the cross section is accounted 
for with the A^* states 5ii(1650), Pii(1710), and Pi3(1720). To explain the broad 
the energy dependence of the cross section seen by SAPHIR 18] and CLAS {20! 
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iV*(1895) (D13) 


* 


* 






A* (1900) (531) 






* 


* 


A*(1910) (P31) 






* 


* 


K*(892) 


* 


* 


* 


* 




* 


* 


* 




A*(1800) (5oi) 




* 






A* (1810) (Poi) 




* 




* 


E*(1880) (Pii) 








* 



TABLE I: Resonances included in the hadrodynamic Bennhold and Mart model [26(] and the 
Janssen et al. model [25'] highlighted in this work for both the K and K'^YP final states. 

BM model includes a spin-3/2 Di3(1895) resonance that was predicted in the relativized 
quark model of Capstick and Roberts to have a strong coupling to the K'^A. channel, but 
which was not well established from existing pion-production data. In addition this model 
includes t-channel exchange of the vector ir*(892) and pseudovector Ki{121{]) mesons. In 
this model the inclusion of hadronic form factors leads to a breaking of gauge invariance 
which is restored by the inclusion of counter terms following the prescription of Haberzettl 



27|. 



For K'^TP production, the BM model includes the Born terms as well as the A^* resonances 
5ii(1650), Pii(1710), and Pi3(1720), and the A* resonances 53i(1900) and P3i(1910). The 
model also includes ii"*(892) and -ft'i(1270) exchanges. The modeling of hadronic form factors 
for the TP channel is handled as described above for the A channel. The BM model does 
not include any u-channel diagrams for either KY final state. 

rn 

In this work, we also compare our data against model B of Janssen et al. [25i], which 
counterbalances the strength from the Born terms by introducing hyperon resonances in 
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the M-channel, where a destructive interference of the w-channel hyperon resonance terms 
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claim that this is 



with the other background terms occurs. The authors of Ref. 
a plausible way to reduce the Born strength. For the K^K calculations, the included u- 
channel resonances are the A* states S'oi(1800) and Poi(1810). For the K+TP calculations, 
the included w-channel resonances are the A* Poi(1810) and the S* Pii(1880). Ref. 29| 
states that there is very little theoretical guidance on how to select the relevant resonances 
and how to determine realistic values for the associated coupling constants. It is stated 
that the same qualitative destructive interference effect was observed for other w-channel 
resonance choices, and that the introduced resonances should be interpreted more properly 
as "effective" particles that account for a larger set of hyperon resonances participating in 
the process. The s-channel and t-channel resonances included in the JB model are nearly 
the same as in the BM model. Hadronic form factors are included in the model with gauge 
invariance restoration based on the approach by Gross and Riska 30 1. 

Different models have markedly different ingredients and fitted coupling constants. Cer- 
tainly not every available hadrodynamic model is discussed in this work. However, it is 
worth mentioning that analysis of Saghai et al. 3l|, using the same data set employed for 
the BM and JB models, have shown that by tuning the background processes involved in the 
K^K reaction in the form of additional w-channel resonances, the need to include the extra 
Di^ N* state was removed. Another analysis including newer photo- and electroproduction 



data from JLab by Ireland et al. 



32l | has shown some evidence for the need for an additional 



N* state at about 1900 MeV (one or more of ^n, Pu, P13, -D13), however they concluded 
that a more comprehensive data set would be required to make further progress. 

A recent coupled-channels analysis by Sarantsev et al. 33|] of the photoproduction data 
from SAPHIR and CLAS, as well as beam asymmetry data from LEPS for K'^A 23|] and 
data from tt and rj photoproduction, reveals evidence for new baryon resonances in the high 
W mass region. In this analysis, the full set of data can only be satisfactorily fit by including 
a new Pu state at 1840 MeV and two D13 states at 1870 and 2130 MeV. Of course these 
fits have certain ambiguities that can be resolved or better constrained by incorporating 
electroproduction data. 

The CLAS and SAPHIR photoproduction experiments measured only the ax term. The 



more recent CLAS data 
to the SAPHIR data [18|, 



2l| , with higher statistical precision and finer binning compared 



19j, reveal that the strength and centroid of the W structure near 
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1.9 GeV changes with angle, indeed pointing to the possible existence of more than one 



s-channel resonance as suggested by the analysis of Ref. [33|]. The interference structure 
functions, (Ttt and (Jlti which are accessible in the electroproduction data and presented 
in this work, will be useful in further constraining and testing models that include new 
s-channel resonance diagrams in this mass region. In addition, the ctt and ctl structure 
functions will also provide crucial information to constrain the model parameters for the 
resonance and background diagrams. 

B. Regge Models 

In this work we also compare our results to a Reggeon-exchange model from Guidal, Laget, 



and Vanderhaeghen [3J] (referred to here as the GLV model). This calculation includes no 
baryon resonance terms at all, but is instead based only on gauge invariant t-channel K and 
K* Regge trajectory exchange. It therefore provides a complementary basis for studying 
the underlying dynamics of strangeness production. It is important to note that the Regge 
approach has far fewer parameters compared to the hadrodynamic models. These include 
the K and K* form factors, which in the GLV model are assumed to be of a monopole form 
Fk,k* = [1 + Q"^ I K\ with a mass scale Kk^k'=^-^ GeV^ chosen to reproduce the JLab 
Hall C (Tl, (Tt data jlJI- In addition, the model employs values for the coupling constants 
Qkyn and Qk'YN taken from photoproduction studies. 

The model was fit to higher-energy photoproduction data where there is little doubt of the 
dominance of these kaon exchanges, and extrapolated down to JLab energies. An important 
feature of this model is the way gauge invariance is achieved for the kaon t-channel exchange 
by Reggeizing the s-channel nucleon pole contribution in the same manner as the kaon t- 



channel diagram 35|]. This approach has been noted as a possible reason why the Regge 
model, despite not including any s-channel resonances, was able to reproduce the JLab Hall 
C o"l/o"t data at Q"^ > 0.5 GeV^. The stated reason is that due to gauge invariance, the 
t-channel kaon exchange and s-channel nucleon pole terms are inseparable and must be 
treated on the same footing. In the GLV Regge model these terms are Reggeized in the 
same way and multiplied by the same electromagnetic form factor. No counter terms need 



to be introduced to restore gauge invariance as is done in the hadrodynamic approach 34 1. 
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V. EXPERIMENT DESCRIPTION AND DATA ANALYSIS 



The experiment was performed using the electron beam at JLab and the CLAS detector 
in Hall B. An electron beam of 5 nA current was incident upon a 5-cm long liquid- hydrogen 
target, resulting in an average beam-target luminosity of £ ~ 10^"^ cm~^s~^. Data were 
taken with beam energies of 2.567, 4.056, and 4.247 GeV. In this analysis the 4.056 and 
4.247 GeV data have been combined into a single data set, referred to throughout this 
work as the 4 GeV data set. The 2.567 data set has a live-time corrected luminosity of 
about 1.32 fb^^, while that for the 4.056 and 4.247 GeV data sets are about 0.67 fb~^ and 
0.80 fb~^, respectively. The beam was effectively continuous with a 2.004 ns bunch structure. 
The large acceptance of CLAS enabled us to detect the final state electron and kaon over a 
broad range of momentum transfer and invariant energy W as shown in Fig. [3l 




FIG. 3: (Color online) CLAS kinematic coverage in terms of vs. W for p(e, e! K'^yV {Y = A, S*^) 
events for beam energies of 2.567 GeV (left) and 4 GeV (right). 



CLAS is a large-acceptance detector j36| used to detect multi-particle final states from 
reactions initiated by either real photon or electron beams. The central element of the 
detector is a six-coil superconducting toroidal magnet that provides a mostly azimuthal 
magnetic field, with a field-free region surrounding the target. The integrated field strength 
varies from 2 T-m for high momentum tracks at the most forward angles to about 0.5 T-m for 
tracks beyond 90°. The field polarity was set to bend negatively charged particles toward 
the electron beam line. Drift chambers (DC) situated before, within, and outside of the 
magnetic field volume provide charged-particle tracking with a momentum resolution of 1- 
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2% depending upon the polar angle within the six independent sectors of the magnet [STJ. 
To protect the chambers from the charged electromagnetic background emerging from the 
target, a small normal-conducting "mini-torus" magnet was located just outside the target 
region. The integral magnetic field of the mini-torus is about 5% that of the main torus. 



t chambers con- 
3811 ■ scintillators 



The outer detector packages of CLAS that surround the magnet and dri 
sist of large-volume gas Cerenkov counters (CC) for electron identification 
(SC) for triggering and charged particle identification via time-of- flight [39|, and a lead- 
scintillator electromagnetic shower counter (EC) used for electron/pion separation as well 
as neutral particle detection and identification [40']. An open trigger for scattered electrons 
formed from a coincidence of the CC and EC signals within a given sector gave event rates 
of about 2 kHz. The total beam charge was integrated with a Faraday cup to an accuracy 
of better than 1%. 

The offline event reconstruction first identified a viable electron candidate by matching 
a negatively charged track in the DC with hits in the SC, CC, and EC counters. The hits 
in the CC and EC counters were required to be within a fiducial region where the efficiency 
was large and uniform. The track was projected to the target vertex to estimate the event 
start time; the estimate was compared to the phase of the accelerator radio-frequency (RF) 
signal to determine this time to better than 50 ps (cr). In contrast to a straightforward 
subtraction of the electron start time from the time, this use of the highly stable RF 
phase improved the hadronic time-of-fiight measurements by almost a factor of \/2. 

A positively charged kaon candidate was identified as an out-bending track found in the 
DC that spatially matched to a SC hit that projected back to the target. The measured 
time-of-fiight of the track and the fitted path length were used to calculate the velocity of 
the particle. This velocity and the measured momentum were used to calculate the mass 
of each charged hadron. For the data discussed here, the kaon momentum range was be- 
tween 300 MeV (software cut) and ~3 GeV (kinematic limit), with a typical fiight path 
of 4.5 m. The measured mass resolution was primarily due to the reconstructed time-of- 
fiight resolution, which was 190 ps (cr) on average, and also included contributions from 
the 1.5% momentum resolution and 0.5 cm path-length uncertainty. A loose cut around 
the reconstructed kaon mass was used to initially select the kaon candidates (a data fil- 
tering condition), however a large background of positively charged pions and protons still 
remained. A momentum-dependent mass cut was used to select the events for the final 
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analysis as shown in Fig. H] for our filtered data files. 




Momentum (GeV) . Mass (GeV) 

FIG. 4: (Color online) Reconstructed mass for positively charged particles. The left figure shows the 
mass plotted against the measured momentum. The lines show the mass cuts used to identify kaon 
candidates. A logarithmic yield density scale is employed. The right figure shows the reconstructed 
mass. These spectra were made from our kaon-filtered data files. 

Corrections to the electron and kaon momenta were devised to correct for reconstruction 
inaccuracies. These arise due to relative misalignments of the drift chambers in the CLAS 
magnetic field, as well as for uncertainties in the magnetic field map employed during charged 
track reconstructions. These corrections were typically less than 1%. 

Using the four-momenta of the incident electron, scattered electron, and candidate, 
the missing mass, corresponding to the mass of the recoiling hyperon, was calculated. The 
missing-mass distribution contains a background that includes a continuum beneath the 
hyperons that arises due to multi-particle final states where the candidate results from a 
misidentified pion or proton, as well as events from ep elastic scattering (protons misidentified 
as kaons) and events from Ti^n final states (pions misidentified as kaons). The elastic events 
are kinematically correlated and show up clearly in plots of 6'|^ (cm. angle) versus missing 
mass and 6k (lab angle) versus (Figs. [5^ and b, respectively). A cut on the elastic 
band in the 6k versus space removes this contribution with a small loss of hyperon yield 
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that is later accounted for with our Monte Carlo generated acceptance function. The vr+n 
events are removed with a simple missing-mass cut in which the detected kaon candidates 
are assumed to be pions. Typical missing-mass distributions showing clear A(1116) and 
S''(1193) peaks are shown in Fig. [6l 




0.6 0.8 1 1.2 1.4 1.6 0.5 1 1.5 



MM(eK"") (GeV) . (GeV^) 

FIG. 5: (Color online) (a) 0^ versus p{e, e'K^)Y missing mass showing ep elastic events and e7r~^n 
events. The vertical bands correspond to ground state A(1116), S'^(1193), and A(1405)/S'^(1385) 
hyperons. (b) 9k versus for p{e, e'K'^)Y showing the ep elastic events and the cut used to 
remove them. 

The continuum from the multi-particle final states that lies beneath the A and S° hyperon 
peaks in our mass spectra is accounted for by a fitting process in which identified pions and 
protons in our unfiltered data files are assumed to be kaons. Missing-mass distributions are 
generated for each assumption in each of our different bins in Q^, W, cos 6*1^, and $. The 
resulting distributions, along with template shapes for the A and S'' hyperons determined 
from Monte Carlo simulations, are fit to the missing-mass spectra using a maximum-log- 
likelihood method appropriate for the low statistical samples in our four-dimensional bins. 
The template shapes for the hyperons were produced from a simulation that included radia- 
tive processes and was matched to the detector resolution. Typical fits of the missing-mass 
distributions are shown in Fig. [61 The final yields in each kinematic bin were determined 
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by taking the number of counts determined from the fits that fell within a mass window 
around the A (1.095 to 1.165 GeV) and S*^ (1.165 to 2.3 GeV) peaks. Hyperon events in the 
tails of the distributions that fell outside of our mass windows were accounted for by our 
acceptance correction function. After removal of all backgrounds, a total of 1.4 x 10^ K^K 
and 6.7 x 10^ K^TP final state events were obtained across the entire kinematic range for 
the 2.567 GeV data set, while 9.7 x lO'' K+A and 4.7 x 10^ events were obtained for 

the 4 GeV data set. 



No. Of Data Events = 297 
No. Of Fit Events = 297,009 
No, Of A Events = 201,1 16 
No, Of Z° Events = 53,486 



Dota Points ■ 
A Fit 



Total Fit 
□ E'Fit 




= 1 ,0 GeV 
W = 1 .850 GeV 
cos(0 = 0,65 
-105° 



1.1 1.12 1.14 1.16 1.18 1.2 1.22 1.24 1.26 

Missing Mass (GeV) 



1.08 1.1 1.12 1.14 1.16 1.18 1.2 1.22 1.24 1.26 

Missing Mass (GeV) 



FIG. 6: (Color online) Signal and background fits from the 2.567 GeV data for the e' miss- 
ing mass spectrum summed over all kinematics (left) and for a typical cos0|^/$ bin (right) at 
(5^=1.0 GeV^ and VF=1.85 GeV to demonstrate the typical fit quality in our data. 



The data were binned in a four-dimensional space of the independent kinematic variables, 
Q^, W ^ cos6^, and $. Table HTl lists the kinematic bin definitions used in the analysis. The 
data at 2.567 GeV consisted of data sets taken with two settings of the main CLAS torus 
field that were combined together. As mentioned above, our 4 GeV data set consisted of 
data acquired at beam energies of 4.056 and 4.247 GeV at the same torus field setting. When 
combining the data at 4.056 and 4.247 GeV, we evolve the cross sections at 4.247 GeV to the 
bin center of the 4.056 GeV data using our model of the cross section. As the two data sets 
are close in energy and e, there is little systematic uncertainty involved in this procedure. 

In this analysis the number of $ bins was 8 for the two backward-most bins of 6'^ and 
12 for the four forward-most 6'|^ bins. The larger number of $ bins for forward and central 
6^ increased the reliability of the $-fits in the presence of the forward beam "hole" of the 
spectrometer, an area of depleted acceptance corresponding to tracks with small laboratory 
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0.8 - 1.3 
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2.3 - 2.8 


W 


8 


1.60 - 1.70 


5 


1.60 


- 1.70 


8 


1.6 - 1.7 


(GeV) 




1.70 - 1.75 




1.70 


- 1.80 




1.7 - 1.8 






1.75 - 1.80 




1.80 


- 1.90 




1.8 - 1.9 






1.80 - 1.85 




1.90 


-2.00 




1.9 - 2.0 






1.85 - 1.90 




2.00 


- 2.10 




2.0 - 2.1 






1.90 - 1.95 










2.1 - 2.2 






1.95 - 2.00 










2.2 - 2.3 






2.00 - 2.10 










2.3 - 2.4 


cos 6*^ 


6 


-0.8 - -0.4 








6 


-0.8 - -0.4 






-0.4 - -0.1 










-0.4 - -0.1 






-0.1 - 0.2 










-0.1 - 0.2 






0.2 - 0.5 










0.2 - 0.5 






0.5 - 0.8 










0.5 - 0.8 






0.8-1.0 










0.8 - 1.0 




8 


cose*j^ < -0.1 








8 


cos 6*^ < 0.2 




12 


cos6'|^ > -0.1 








12 


cos6'|f > 0.2 



TABLE II: The number of bins and the bin sizes employed for the 2.567 and 4 GeV analyses in 
this work. For the 2.567 GeV data analysis two different binning choices were made for W for each 
bin in Q^, cosO"^, and 

angles. Bins significantly overlapping this forward hole were excluded from our analysis. 

The average differential cross section for each hyperon final state in each bin i was com- 
puted using the form: 
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dai I f 1 \ f Ni \ f 1 



where Ni is the hyperon yield, Ai is the acceptance, A^^o is the hve-time corrected incident 
electron flux, Ri is the radiative correction factor, A^^ is Avogadro's number, p is the target 
density ((p) = 0.072 g/cm^), t is the target length, and A^ is the atomic weight of hydrogen 
(1.00794 g/mol). The product AQ'^AWAcos9*j^A^ represents the volume of the z*^ bin 
corrected for kinematic limits. 

The geometric acceptance and reconstruction efficiencies were calculated using a standard 



model of the CLAS detector based upon a GEANT simulation [4ll|. To reduce the model 
dependence of the computed CLAS acceptance, it is important to match the distributions of 
accepted data and Monte Carlo events as a function of the relevant kinematic variables Q^, 
W, cos^l^, and $ for both the K^A and final states. This match must be ensured at 

all beam energies and torus field settings employed in the analysis. 

A variety of reaction models (see e.g. Ref. were employed as input event gener- 
ators for the simulated events. Because none of the models agreed particularly well with 
our K~^A or K~^Ti^ data, we developed our own models that were able to match the data 
reasonably well over our full kinematic phase space. We developed two different ad hoc 
models for the K^A analysis that both represented our data equally well, although with 
slightly different dependencies on Q^, W, and cos 9}^. One ad hoc model was developed for 
the i^+S^ analysis. These models were used as input to determine our detector acceptance 
function, radiative corrections, and bin-centering correction factors. In addition we devel- 
oped an event generator based on fits to our K~^A data. Differences between our two ad hoc 
models and our data-fitted model for K^A were used to estimate the model dependence of 
our results. Further details are included in Section IVIIi Figure [7] shows the dependence of 
the K^A acceptance upon 6^ and $ for a bin in and W. Typical acceptances of CLAS 
for the e'K^ final state were in the range of 1 to 30% depending on kinematics. Note the 
strong variation in acceptance as a function of both cos 6"^ and $ due to the geometry of 
CLAS. 

In the detector simulation, particles generated at the target were propagated through the 
CLAS magnetic field and were permitted to interact with materials and to undergo decay. 
These tracks then generated simulated detector hits. Hits corresponding to the known dead 
areas in the detectors were removed, and the hits were smeared according to the known 
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FIG. 7: Distribution of the computed K^A acceptance for CLAS as a function of cos 0^ and ^ for 
the W=1.85 GeV and Q^=0.65 GeV^ bin. The depleted area near ^ of 0° for forward angles is the 
forward "hole" in CLAS due to the beam pipe. The number of bins in ^ is different for forward 
and backward cos 0^ . The statistical error bars from the Monte Carlo are smaller than the symbol 
size on this plot. The curves on each plot serve only to guide the eye. 



detector resolution effects. These simulated events were passed through the same analysis 
chain as the real data. Geometrical fiducial cuts were applied to both the data and simulated 
events to eliminate areas of inefficient detector response or where the response was not well 
modeled. These areas were typically within a few degrees of the magnet coils and near the 
edges of the Cerenkov detector. 

The radiative correction for each kinematic bin was computed from the ratio of the model 
cross sections with and without radiative effects. We used two very different methods to 
compute this correction factor. 

The first method used an acceptance-rejection technique, where events were generated 
uniformly in W, Q^, cos 6*^, and $, with a weight determined via the cross section model of 
Ref. Q]. The energies of externally radiated photons from the incident electron, and from 
the emitted electron and kaon in the region of the target proton, were generated according to 



the formulas of Mo and Tsai 



42]. The weight of the event was adjusted to account for hard 



and soft internal radiative effects, and the post-radiation kinematic variables were calculated 
to identify the bin into which the event fell. While lacking in computational efficiency, this 
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method benefited from being able to compute cross sections as well as the expected event 
distributions, both with and without radiative effects, thus providing a consistent event 
sample to use for the acceptance studies. The second method used the same formula for 
calculating the radiated cross section from the assumed non-radiated model. However, it 
integrated the resulting six-dimensional cross section over the two unseen dimensions that 
corresponded to a radiated photon from either the initial or scattered electron. The ratio 
of the integrated radiated cross section (now corresponding to a four- dimensional space) to 
the unradiated cross section yielded the radiative correction factor. We chose to use this 
second method because of its superior computational speed. It was extensively checked 
versus the "EXCLURAD" code of Afanasev [4^ . Differences between the two methods allowed 
us to estimate the size of any residual uncertainty in the radiative correction procedure (see 
Section Ivni). 

In order to do a full separation into four structure functions, we can fit our full set of data 
including the differential cross sections from both beam energies with a function of the form 
/($, e); the fitted parameters being the values of 0"^, cl, cttt, and 0"^^ at some fixed point in 
Q^, W, and cos^|^. Alternatively, we can extract au in each Q^, W, and cos^|^ bin from a 
$ fit for each beam energy separately, and then do a linear e fit to separately extract ctt and 
(Jl; this being the well-known Rosenbluth separation technique. In either case, in order to do 
these fits, we must first define the cross sections at a specific fixed point within the bin, and 
not merely as an average over a given bin volume. This is especially true when the bins are 
large and event-weighted average values of kinematic variables can be different for different 
$ bins. Using an integration over our model cross section, we calculate the cross section 
at the fixed point given its average over the bin volume. This correction is referred to as 
a "bin-centering correction", or more accurately, as a "finite bin size correction". By using 
our model event generators, we simply calculated the ratio of the cross section evaluated at 
the assigned bin center to the average cross section integrated over the bin to obtain the 
finite bin size correction factor. Systematic uncertainties associated with these corrections 
were extensively studied (see Section Ivnll . 

In Fig. [8] we show a sample of the ^-dependent differential cross sections for the K~^\ 
final state at representative kinematic points. The different shapes of the differential cross 
sections vs. $ in each of our bins in Q^, W, and cos^^ reflect differences of the interference 
terms, (Ttt and ctlt- The differences in scale reflect the differences in au. 
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FIG. 8: <I>-dependent differential cross sections (nb/sr) and fits for K events from our 2.567 GeV 
data at ^^=0.65 GeV^ for eacfi of our six cos 0^ bins (labeled at the top of each column) for three 
different W bins (labeled on the left of each row). The curves represent fits to the ^-dependent 
differential cross sections. 

VI. STRUCTURE FUNCTION EXTRACTION 



The full set of differential cross sections da-^/dTL*^ included in this work for each hyperon 
final state consists of 156 bins in Q^, W, and cos^l^^ for the 2.567 GeV data (accounting for 
both W binning scenarios given in Table [Tll) and 192 bins for the 4 GeV data. This amounts 
to 1664 data points in Q^, W, cos6^, and $ for the 2.567 GeV data and 1920 data points at 
4 GeV for each hyperon final state. In this section we provide details regarding the structure 
function extraction. In Section IVI Al we focus on the separation of the structure functions 
au, (Jtt, and air- In Section IVIBI we present the extraction of ax and cr^ separately from 
(7(7 using a Rosenbluth fit and also a simultaneous e — $ fit of our data at 2.567 and 4 GeV. 
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A. Extraction of ajj, (Jtti and a^j^ 



The differential cross sections plotted in Fig. [8] are actually the mean values within the 
finite size of the $ bins and therefore do not necessarily reflect the value at the bin center. 
Thus directly fitting these data with eq.© to extract the structure functions au = (JT + eaL, 
aTT, and clt would be inappropriate. Integrating eq.© over the finite bin size, A$ = 

— where and $i are the upper and lower limits of the bin, respectively, gives: 



(To = -r-r- / ( cr;7 + ecTTT COS 2$ + A/e(e + 1)<7lt cos $ j (i$ 

= ^ (^(T[7 A$ + I^TT (sin 2$„ - sin 2$^) + ^ e{e + 1)(Tlt (sin - sin j . (7) 

(To now represents the value of the measured bin-averaged cross section in a given $ bin and 
fitting the data with eq.([7]) yields the separated structure functions for a given bin in Q^, 
W, and cos6'^^. The "e" pre-factors were evaluated at the bin center and divided out. 

Prior to the $ fits, the statistical uncertainty of each cross section bin was combined 
linearly with that portion of the systematic uncertainty arising from the yield extraction 
procedures (see Section [VIII for details). A few points were removed from the fits based 
upon their low acceptance in CLAS, in order to prevent bins with a very small acceptance 
from distorting the extracted structure functions. A point was rejected if its acceptance at 
2.567 GeV (4 GeV) was less than 2.0% (1.0%) or less than 10.0% (5.0%) of the average 
acceptance over all bins at the same Q^, W, and cos6'^. 

In reporting the final results from our $ fits, several cos^l^^ bins have been discarded. 
In general, these bins were near the edge of our kinematic acceptance and had limited $ 
coverage. In addition, the statistical uncertainties were large on the points in these bins that 
survived the acceptance criteria described above. Typically, the missing points were near 
$ = ±7r or $ = - exactly where points are needed to constrain the interference structure 
functions. The resulting $ fits for these bins had x^/^ values, where u represents the number 
of degrees of freedom, that were uniformly too small considering the expected x^/ ^ distri- 
butions. In other words, a three-parameter fit of these bins had too many parameters, given 
the low number of data points and the large uncertainties, to give unambiguous solutions 
for the structure functions. We also examined the x'^/'^ distributions for the remaining fits 
and found that they were well represented by their expected probability distributions, which 
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instills confidence in the quality of the data, the assigned uncertainties, and the fits. 
B. Separation of ax and 

The extraction procedure detailed in Section IVI Al yielded the bin-centered structure 
functions au, cttt, and air- To further separate au into its component parts, ut and o"l, 
we have two options. The first is the standard Rosenbluth separation technique, in which 
au is determined for two different beam energies (or different e values) but for the same 
point in Q^, W, and cos 6*^, and fit as a linear function of e. An alternative approach is to 
simultaneously fit the data from the two energies as a function of e and this time explicitly 
replacing au in eq.Q with ar + €aL- This method has the advantage of constraining the 
individual parameters, ax, a^, axr, and air to have the same value for the two different 
beam energies, as they must since they are explicit functions of Q^, W, and 6^ only. This 
approach represents an important systematic check as the forward beam hole of CLAS affects 
the acceptance function differently at 2.567 GeV relative to 4.056 and 4.247 GeV. 

The separation of the structure functions ar and ai can only be performed in {Q^, W, 
cos 6*1^) bins where the 2.567 GeV and 4 GeV data overlap. Figure [9] shows plots of e vs. 
for four different 100-MeV wide bins in W from 1.65 GeV to 1.95 GeV, and highlights the 
kinematic coverage of CLAS. The cut-off at low is due to the minimum 6e detectable by 
CLAS and the low e cut-off is due to the maximum detectable by the Cerenkov detectors. 
For this analysis, the data overlap only for a rather narrow region at about 1 GeV^. 
We have performed a separation of ax and ai for Q'^=1.0 GeV^ for the i^+A final state for 
W=1.65, 1.75, 1.85, and 1.95 GeV, and at values of W=1.75, 1.85, and 1.95 GeV for the 
i^+S^ final state. In this two beam energy separation, we have typical differences in e of 
about 0.4. Of central importance in this analysis is the fact that this separation is performed 
for the first time away from the condition of parallel kinematics (i.e. 61^=0° or along the 
virtual photon direction). 

Before the separation of au could proceed, we first had to account for the binning dif- 
ferences between the 2.567 GeV and 4 GeV data sets. Due to consideration of statistics 
in the two separate data sets, the 2.567 GeV data were sorted in 50-MeV wide W bins for 
the extraction of au, axr, and a^T, while the data at 4 GeV were sorted in 100-MeV wide 
W bins (see Table HI])- In order to perform either the Rosenbluth fit or the simultaneous 
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FIG. 9: CLAS acceptance of the scattered electron in terms of e vs. at 2.567 GeV (lower bands) 
and 4 GeV (upper bands) for four 100-MeV W bins centered from 1.65 to 1.95 GeV. The vertical 
line at Q^=1.0 GeV^ marks where we have performed the separation of ax and a^. 



e — $ fit, the 2.567 GeV data had to be resorted into W bins that were 100-MeV wide. 
In computing the cross sections for the 100-MeV wide W bins at 2.567 GeV, the hyperon 
yield fits were redone and all other factors associated with computing the cross section were 
recalculated using Monte Carlo based on the 100-MeV wide bin. 

The Rosenbluth extraction procedure is a standard technique to separate ax and 0"^. The 
error bars on these structure functions result from the statistical and systematic uncertainties 
on the two au cross section points used in the extraction. With only two data points, 
the slope parameter (ctj^) and the intercept parameter (ctt) along with their associated 
uncertainties, can be computed analytically. Figure [10] shows a representative plot of the 
au cross sections for the K^A final state at W=1.85 GeV for each of our six cos6'|^ bins. 
This plot also serves to indicate the typical e values and spread for the two data sets. The 
data points at 2.567 GeV and 4 GeV have each been evolved using our ad hoc models to 
the point Q^=1.0 GeV^. The analysis employs the highest bin from our 2.567 GeV data 
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set and the lowest bin from our 4 GeV data set. 
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FIG. 10: Representative Rosenbluth separation plots of au (nb/sr) vs. e for our K^K data at 
(5^=1.0 GeV^ and VF=1.85 GeV for our six cos6*^ bins. The lines represent the fits that determine 
the slope parameter {(Jl) and the intercept parameter (ctt), which are printed on each plot. The 
error bars on the data points and the errors listed for ai and ctt on each plot represent the combined 
statistical and systematic uncertainties. 



An example of the comparison between the separate $ fits for the 2.567 GeV and 4 GeV 
data and the simultaneous fit for both energies for the K^K and K^TP reactions is shown 
in Fig. [TTlfor Q^=1.0 GeV^ and 1^=1.85 GeV. The differences between the differential cross 
sections for the 2.567 and 4 GeV data for a given bin (see Fig. [11]) in cos6'|f are due not 
only to the beam energy dependent e pre-factors (defined in eq.([5])), but also to the different 
systematic variations associated with the acceptance functions of CLAS at these energies. 
Of importance is that the simultaneous fits differ from the single beam energy fits only 
where the single beam energy fits have large error bars (e.g. 4 GeV back-angle K^TP bin) 
or are missing $ points due to our minimum acceptance cut-off criteria (e.g. 4 GeV back- 
angle K^K bin). In these cases the simultaneous fit procedure leads to extracted structure 
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function with reduced uncertainties compared to the single beam energy fits. 
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FIG. 11: Comparison of $ fits to the differential cross sections performed for two different al- 
gorithms. In the first approach (dashed curves) the cross sections at 2.567 and 4 GeV are fit 
separately. In the second approach (solid curves) the two different beam energy data sets are fit 
simultaneously. These plots are for kinematics with Q^=1.0 GeV^ and W=1.85 GeV for the K~^A 
(left columns) and iiT+E^ (right columns) final states at different values of cos 6^ as indicated. 



VII. SYSTEMATIC UNCERTAINTIES 
A. Overview 

To obtain a virtual photoabsorption cross section, we reconstruct events with an outgoing 
electron and K'^, and then fit the missing- mass spectra for each of our bins in Q^, W, cos 9"^, 
and $ to obtain the yields for the reactions K'^A and K'^TP. The yields are corrected for 

the acceptance function of CLAS, radiative corrections, and finite bin size effects. Finally, 
we divide by the virtual photon flux factor at the bin center, the bin volume corrected 
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for kinematic limits, and the beam-target luminosity to yield the cross section. Each of 
these procedures is subject to systematic uncertainty. We typically estimate the size of 
systematic uncertainties by repeating a procedure in a slightly different way (e.g. by varying 
a cut parameter within reasonable limits or by employing a slightly different algorithm) and 
noting how the results change. The difference in the results is then used as a measure of the 
systematic uncertainty. In this section we describe our main sources of systematics. 

With respect to their effect on our results, there are three types of systematic effects: 
uncertainties that affect the yield extraction in a seemingly random fashion where the sys- 
tematic uncertainty is proportional to the size of the statistical uncertainty, "scaling" uncer- 
tainties that affect both the cross sections and structure functions by a simple scale factor, 
and ^-dependent uncertainties such as using an event generator with a ^-dependence that 
does not quite match the data. 

These uncertainties are handled in different ways. Because the size of the "yield extrac- 
tion" uncertainty depends on the size of the statistical uncertainty, we take the random-type 
systematic uncertainties into account by enlarging the statistical uncertainty before the $ 
fit to extract the structure functions. These fractional systematic scaling uncertainties are 
multiplied by the value of the cross section or the structure function in question to get the 
absolute uncertainty. The remaining uncertainties, which can in general have $ or other 
kinematic dependencies, are estimated by extracting the ^-dependent structure functions 
for two, similar procedures. This method gives an absolute estimate for a structure function 
uncertainty. 

The primary sources of systematic uncertainty for this experiment came from the Monte 
Carlo model dependence (acceptance, radiative corrections, finite bin size corrections), de- 
tector efficiency, and the yield extraction. With the very large acceptance and a four- 
dimensional kinematic space, systematic uncertainties were studied on a bin-by-bin basis. 
Table HIT] summarizes our estimates of the average systematic uncertainties on the differential 
cross sections associated with various effects. The different types of systematic uncertainties 
mentioned above are referred to as "stat.", "scaling", and ""I'-dep." in the column labeled 
"Type" in Table [ml 

The main categories of systematic uncertainty in this analysis include: (i). event recon- 
struction efficiency (SaER), (ii). yield extraction (Sctye), (hi), model dependence (SaMon), 
(iv). radiative correction theory uncertainty (^ct/j^), (v). virtual photon flux {6a fi^^), and 
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Category 


Type 


Sources 


Avg. Size 


i). 


Event Reconstruction 


scaling 


Trigger+tracking efficiency 


1% 






$-dep. 


Electron fiducial cut 


3.6% 






$-dep. 


Kaon fiducial cut 


4.1% 






scaling 


Electron PID efficiency 


L5% 






scaling 


Kaon PID efficiency 


L0% 






scaling 


CC efficiency 


2-5% 






scaling 


CLAS forward angle response 


1-10% 


ii). 


Yield Extraction 


stat. 


Signal templates 
Background removal 


25% • Stat 


iii) 


. Model Dependence 


<^-dep. 


Acceptance calculations, 
Radiative corrections, & 
Finite bin size corrections 


8.0% 


iv) 


RadCorr: Theory 


scaling 


Integration vs. EXCLURAD 


3.4% 


v). 


Photon Flux-factor 


scaling 


Momentum and angle uncertainties 


3.0% 


vi) 


Luminosity 


scaling 


Live time correction 


0.5% 






scaling 


Faraday cup accuracy 


1.0% 






scaling 


Hydrogen target thickness 


3.0% 



TABLE III: Sources, types, and average sizes of systematic uncertainties on the differential cross 
sections. 

(vi). luminosity {6ac)- Each of these categories is explained in more detail in the next 
section. The final systematic uncertainty assignment to our extracted structure functions is 
explained fully in Section IVII CI While the yield extraction systematic uncertainty, as ex- 
plained below, is treated as an effective increase in our statistical uncertainty, the remaining 
systematic sources are added in quadrature to arrive at our final uncertainty assignment as: 

1 /2 

+ 54c + -^4- + -^4) • (8) 
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B. Systematic Uncertainty Categories 



(i) . Event reconstruction efficiency: This efficiency is a convolution of the charged particle 
track reconstruction efficiency in CLAS, the efficiency of our particle identification algorithms 
for the electron and kaon, and the triggering efficiency. The CLAS trigger and tracking 
efficiency (which are essentially 100%) have been studied and represent small contributions 
to our systematics. The definitions of the electron and kaon fiducial cut boundaries (which 
cut ~10% of our event sample) and the particle identification (PID) cuts (which cut ~15% 
of our event sample) have been varied within reasonable limits to determine their effect on 
the resulting cross sections. Each of these systematic sources is relatively small, and overall 
they contribute about 6% to our total systematic uncertainty. Each source is independent 
of the kinematics of the final state particles. 

There are two additional sources of systematic uncertainty in this category that have a 
value that depends on the final state kinematics. One of these sources accounts for unphysical 
small-scale fluctuations in the measured efficiency function of the Cerenkov detector (which 
has typical efficiencies of 95% in our fiducial region), which were much more apparent at 
forward angles in CLAS. This "CC efficiency" systematic has been assigned as 5% for the 
lowest Q'^ bin for each beam energy data set (Q^=0.65 GeV^ at 2.567 GeV and Q'^=1.0 GeV^ 
at 4 GeV) where the electrons populate smaller angles in CLAS. For all other bins the 
systematic has been assigned to be 2%. The other kinematics-dependent systematic arises 
due to the fact that our ai/ax extraction was performed in a region with only modest 
kinematic overlap between the 2.567 GeV and 4 GeV data sets, namely Q^=1.0 GeV^ (see 
Fig. [9]). The electrons in the 2.567 GeV data sample populate a well understood and well 
modeled portion of the CLAS detector. However the electron sample in the 4 GeV data 
populate the forward-most portion of CLAS where the acceptance is difficult to model due 
to the forward beam hole of CLAS and where the Cerenkov efficiency varies rapidly due 
to the mirror geometry of the detector 38j. With the 4 GeV data at Q^=1.0 GeV^, we 
have assigned a W^-dependent systematic uncertainty that is 10% at W=1.65 GeV, 5% at 
W=1.75 GeV, 2% at W=1.85 GeV, and 1% at 1^=1.95 GeV. 

(ii) . Yield extraction: As discussed in Section |Vl we use Monte Carlo templates that 
have been matched to the data for the A and S° peaks and background forms based on the 
spectra of misidentified pions and protons in order to fit the hyperon missing-mass spectra. 



34 



We studied various changes to our procedures such as changing the histogram bin size in 
the fitting procedure and using different forms for the background shape (e.g. using both 
misidentified pions and protons, only misidentified pions, and only misidentified protons) 
and conclude that all systematic effects get larger in direct proportion to the size of the 
statistical uncertainty. We estimated that any remaining systematic uncertainty due to the 
yield extraction is roughly equal to 25% of the size of the statistical uncertainty in any given 
bin. We added these correlated uncertainties linearly with the statistical uncertainties on 
our differential cross sections before performing the $ fits. 

(iii). Model dependence: We have studied the systematic uncertainty associated with 
the model dependence of the convolution of the CLAS acceptance correction, the radiative 
corrections, and the finite bin size correction together because they are correlated, especially 
by their sensitivity to the underlying physics model that we use for the Monte Carlo event 
generator. Specifically we studied the overall model dependence by varying the physics 
model used in our Monte Carlo program and stepping through the full analysis chain from 
yields to cross sections to structure function extraction. 

We tried a number of existing hadrodynamic models, but found the agreement with our 
data to be unsatisfactory. Ultimately we employed the model of Bennhold and Mart j26| and 
adjusted the parameters in an ad hoc fashion to get a better match to our measured K^h. 
and K'^TP cross sections as a function of Q^, W , cos9'^, and $ (see discussion in Section IVT) . 
In Fig. [12] we show comparisons of our ad hoc event generator models to our initial model 
from Bennhold and Mart 26|] and to our data. Our studies showed that the event-generator 



model dependence introduced an average systematic uncertainty on our differential cross 
sections of 8%. 

(iv) . Radiative correction theoretical uncertainty: The radiative correction factor was 
calculated using a multi-dimensional integral approach (see Section |V]). To calculate the 
theoretical uncertainty, our results were compared to the exact one-loop calculations from 
the EXCLURAD code [44]. The average deviation was approximately 3.4% over all kinematical 
bins. 

(v) . Virtual photon flux factor: We estimated uncertainties on the average virtual photon 
flux factor across our kinematics by propagating through the flux definition (see eq.(l3])) 
the uncertainties associated with W and that arise from the absolute uncertainty in 
the reconstructed electron momentum and angles. The uncertainty in the fiux factor was 
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FIG. 12: Cross section models for the K'^ K (left) and K'^TP (right) structure functions (in units 
of nb/sr) vs. W compared to our CLAS data (square points) at cos6'^ = 0.05 and 0.35. The ad 
hoc models employed (dashed and dot-dashed curves - discussed in Section |V]) were based on the 
Bennhold-Mart model (solid curve) as a starting point. The plots are for the 2.567 GeV data 
set at Q2=o.65 GeV^. 



determined to be less than 3%. 

(vi). Luminosity: The uncertainty in our luminosity is based on the uncertainty in our 
electron flux, target thickness, and measured live time. The total systematic uncertainty 
from these sources is assigned as 3.2%. 



C. Final Systematic Uncertainty Assignments 

The relative systematic uncertainties on the interference structure functions 5(Jtt/c"tt 
and daLx/o'LT must be interpreted with some caution as both of these interference struc- 
ture functions are frequently small in our kinematics. In this regard, defining a relative 
uncertainty is mathematically meaningless. We have chosen instead to quote all systematic 
uncertainties relative to au- Fig. [13] shows that the kinematic-independent systematic un- 
certainties on each of the structure functions au, cttt, and air relative to au are reasonably 
independent of Q^, W, cos 9'^, and $. For this reason we have decided to quote the relative 
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systematic uncertainty as the mean of these distributions for each beam energy. This ehm- 
inates the fluctuations in the determination of the systematic uncertainties associated with 
low statistics portions of our phase space. From these distributions we compute the mean 
and then add in quadrature the systematics associated with the Cerenkov detector efficiency 
(a Q^-dependent systematic) and the forward angle response of CLAS (a l^-dependent sys- 
tematic) to get the final total systematic uncertainty assigned to our data points. The 
same systematics determined from the analysis of the K^K final state are assigned to the 
data for the K^TP final state as the K^K data has smaller statistical uncertainties. The 
final total systematic uncertainty assignments relative to ctu for our three structure function 
separations are given in Table IIV[ 







Systematic Uncertainty 


Beam Energy 


Term 


g2=0.65 GeV2 


Q2=1.00 GeV^ 


9 KfiV Pol/ 

z.oD/ Lxe V 


bau jau 


9.6% 


8.4% 




5aTT / <^u 


11.7% 


10.8% 




Sctlt/ctu 


7.8% 


6.3% 


Beam Energy 


Term 


Q2=i.oo GeV2 


Q2=i.55, 2.05, 2.55 GeV^ 


4 GeV 


6au/cru 


W=1.65 GeV 
1.75 GeV 
> 1.85 GeV 


13.9% 
10.8% 
10% 


8.4% 






W=1.65 GeV 
1.75 GeV 
> 1.85 GeV 


15.4% 
12.7% 

12% 


10.8% 




SaiT/cru 


W=l.65 GeV 
1.75 GeV 
> 1.85 GeV 


12.7% 
9.3% 

8% 


6.3% 



TABLE IV: Total systematic uncertainties assigned to our structure function measurements for 
both the K^A and AT^S*^ final states as a function of kinematics. Uncertainties for au, (Ttt, and 
are all quoted relative to ajj. 

The systematic uncertainty analysis on the separated structure functions cr^ and was 
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FIG. 13: (Color online) Total kinematic-independent systematic uncertainties from eq.dS]) on the 
structure functions au, aTT, and aLT normalized to au for the K^A data as a function of bin 
number for the 2.567 GeV (top) and 4 GeV (bottom) data sets. The wide vertical boundaries 
indicate the bins, the narrow vertical boundaries indicate the W bins within each range, 
and the six points in each W bin represent the angle bins from cos0^=-O.65 to 0.90. 



carried out only for the Rosenbluth separation method. In order to be conservative, the same 
systematic uncertainty is assigned to ctt and ai extracted from the simultaneous e — $ fit. 
This was done as we could not fully disentangle the point-to-point and scale-type systematic 
uncertainties between the two beam energy data sets in the e — $ fits. In this analysis, we 
simply use the two different techniques as a way to perform a consistency check on our 
extracted structure functions. Our analysis shows very good agreement between the two 
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techniques giving us confidence in our assigned systematics. 

We performed several consistency cliecks on our data. Tlie most important were tliat 
cross sections at 2.567 GeV taken witli two different magnetic field settings and our cross 
sections taken at 4.056 and 4.247 GeV agreed within the quoted systematics. This tested 
the accuracy of our knowledge of the acceptance because it varied strongly with field setting 
and beam energy. The other check was to fit the two beam energy data sets simultaneously 
in each of our bins to verify that the relative normalization factor between the two data sets 
was consistent with unity. 



VIII. RESULTS 



A. Three Structure Function Separation 



1. Angular Dependence 

In Figs. [TH and [15] we show the extracted structure functions au, cttt, and air versus 
cos^^ for K+A and K+T.'^ for different W points at Q2=o.65 GeV^ from our 2.567 GeV data 
set. Although we focus on the cos 6*1^ dependence of our low data set at 2.567 GeV, the 
general conclusions that can be drawn from studying the angular dependence are similar for 
our other data sets. However, the full set of our data is available in Ref. 51]. In these plots, 
the data are sorted into W bins 50-MeV wide, except for the first and last W bins which 
are 100- MeV wide (see Table HTll . All data points have been evolved to the given Q^, W, 
and cos 6}^ bin centers. The curves shown are from the hadrodynamic models of Bennhold 



and Mart (BM) [45| (dot-dashed curves) and of Janssen et al. (JB) [46] (solid curves), and 
the Reggeon-exchange model of Guidal et al. (GLV) [47] (dashed curves). 

A number of observations can be made independent of the model calculations. First, 
we observe that the K~^A and K~^Ti^ electroproduction dynamics are very different. The 
data in Figs. [H] and [15] reveal that au is more forward-peaked in i^+A production than 
for production across our full range in W. With regard to the interference structure 

functions, cttt for K^A is roughly one-fourth of au and always negative and very similar 
in structure and magnitude to ctlt, while cttt for K'^'E^ is generally smaller in magnitude 
than for K^A with a peaking at more mid-range angles. The K^A reaction has a significant 
aiT component in the forward direction compared to au, while air for the K^T,^ reaction 
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FIG. 14: Structure functions au, aTT, and aiT (in nb/sr) for K^A production vs. cos9'^ at 2.567 GeV for Q^=0.65 GeV^ and W from 
1.650 to 2.050 GeV. The error bars represent the statistical uncertainties only. The relative systematic uncertainties to ajj are given in 
Table IIVI The curves shown are from the model calculations of Bennhold and Mart (BM) 



45l | (dot-dashed), Janssen et al. (JB) [46| (solid), 



and Guidal et al. (GLV) 47[ (dashed). The models are described in the text. 
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FIG. 15: Structure functions au, <7tt, and aLT (in nb/sr) for K~^T,^ production vs. cos^|^ at 2.567 GeV for Q^=0.65 GeV^ and W from 
1.725 to 2.050 GeV. The error bars represent the statistical uncertainties only. The relative systematic uncertainties to ajj are given in 
Table IIVI The curves shown are from the model calculations of Bennhold and Mart (BM) 



45l | (dot-dashed), Janssen et al. (JB) [46| (solid), 



and Guidal et al. (GLV) 47[ (dashed). The models are described in the text. 



is everywhere consistent with zero. 

The forward-peaking of au and ctlt for K^K compared to K^TP can be quahtatively 
explained by the effect of the longitudinal coupling of the virtual photons. We note that the 



two channels are of nearly equal strength at Q^=0 GeV^ 20|,l2l|, while here at (^^=0.65 GeV^ 
the K^K channel is stronger than the K^TP channel at forward angles by a factor of 2 to 
3. For transverse (real) photons, the t-channel mechanism at low t is dominated by vector 
K*^ exchange, which relates directly to the relative magnitudes of the gK*YN coupling 
constants to the Qkyn constants. As rises from zero, the photon can acquire a longitudinal 
polarization and the importance of pseudoscalar exchange increases. Given that g\j^^ 3> 
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49j, this effect increases the cross section for K^K relative to K'^TP . This 



argument was already noted in the earliest reports of hyperon electroproduction [8|, and 
is strengthened by our observation of a sizeable (Tlt for K^K and a air consistent with 
zero for K'^TP. It should also be the case that since gK*j:N ^ gKT.Ni K* exchange should 
dominate the S° channel. Because K* exchange must vanish at forward angles due to angular 



momentum conservation, the S*^ cross section should also decrease at forward angles 34l |. 

None of the three different models shown is particularly successful at describing all of 
the data. In general the models better agree with the K~^A data than with the data. 
The three models tend to reproduce the qualitative fall-off in cos^^ of au for the K~^A 
data but do not include sufficient forward-angle strength for W < 1.8 GeV. At higher W 
the BM model generally reproduces au, while the JB model consistently is too large at 
forward and backward kaon angles. The GLV model goes above our data as cos 6"^ ^ 1, but 
describes the structure of the K'^A data surprisingly well considering that it has no built-in 
s-channel resonances, au for the K^T.^ data is poorly described by all models, especially 
the JB model, which includes too much w-channel strength, while the BM and GLV models 
generally include too little strength or miss the broad peaking about cos 6^ ~ 0. 

Within the GLV Regge model, the functions axr and auT arise from the interference of 
the K and K* Regge trajectories. This modeling is sufficient to qualitatively reproduce the 
behavior of both the K^A and K^T.'^ data over our full kinematic phase space. The quality 
of the comparisons of the hadrodynamic models to the cttt and auT data are much less 
favorable. For the JB model, axr for K~^A has the correct sign, but its strength is too small 
and the angular dependence does not match the data. For K~^T.^, the JB model predicts 
arprp ~ everywhere, in strong disagreement with the data. For the BM model, axr for K^A 
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has a strength and angle dependence that quahtatively matches the data, but has the wrong 
sign. For K^TP, the BM model has the wrong sign for (Ttt and doesn't match the angular 
distribution of the data. From the JB model, a^T for K^K is consistent with zero at low W , 
but increases in strength for higher W , where the model has the wrong sign compared to 
the data. For K'^TP , the JB model has both the wrong sign and angular dependence. For 
the BM model, (Tlt follows the trends of the K^K data but has overall too little strength, 
while it is reasonably consistent with the K^TP data. 



2. Energy Dependence 

Even if A production for forward-going mesons is dominated by t-channel exchange, 
there is still room for s-channel resonance contributions at more central angles and at all 
angles for the TP . To more directly look for s-channel resonance evidence, the extracted 
structure functions are presented as a function of the center-of-mass energy W for our six 
bins in cos 6^1^. Figures [T6] and [TTI show the results for our 2.567 GeV data at (5^=0.65 GeV^ 
for the contiguous angle bins centered at cos6'^ = -0.6, -0.25, 0.05, 0.35, 0.65, and 0.90. 

Several characteristics of the data stand out. For K^K production, au shows a broad 
peak at about 1.7 GeV at forward angles, and two peaks separated by a dip at about 
1.75 GeV for our two backward angle points. Across our phase space, cttt and aiT for 
K^K production are predominantly negative and about one-third the size of au- Where 
the statistical uncertainties on our data are reasonable (away from the most forward-angle 
point), (Ttt and air seem to be similar in shape to cT[/, but opposite in sign. The K^TP 
structure functions have a different set of features. Both au and axr exhibit a broad bump 
at about 1.85 GeV, while cttt is consistent with zero everywhere. The TP shapes are similar 
for both forward and back-angle production, with a strong peaking at central angles. 

We argue here that our spectra likely reflect the existence of a few underlying s-channel 
resonances along with t-channel processes, but acknowledge that the physical interpretation 
is not straightforward and will require detailed modeling. The ly-dependence of the K^A. 
data for (Ju near threshold shows more structure than a model based upon only t-channel 
exchanges provides (GLV model - dashed curves) and is probably evidence of resonance 



activity. In this range of W , the S'ii(1650) is believed to be dominant in the s-channel 26|. 



There are also a number of known A^* resonances near 1.7 GeV that can contribute to the 



43 



400 
350 



COS ei=-o.6o 



cos 9i=-0.25 



cos 6i=0.05 



cos 



cos 9i=0.65 



cos 61=0.90 





300 T 








250 i- 








200 - 








150 - 


b 






100 r 




1.6 1.7 1.8 1.9 2.0 2.1 1.7 1.8 1.9 2.0 2.1 1.7 1.8 1.9 2.0 2.1 1.7 1.8 1.9 2.0 2.1 1.7L 1.8 1.9 2.0 2.1 1.7 1.8 1.9 2.0 

W(GeV) W(GeV) W (GeV) W (GeV) W (GeV) W (GeV) 



2.1 



FIG. 16: Structure functions ajj, <Jtt, and aiT (in nb/sr) for K^A production vs. W at 2.567 GeV for (5^=0.65 GeV^ for our 6 cosO'^ 
bins. The error bars represent the statistical uncertainties only. The relative systematic uncertainties to au are given in Table llVi The 
curves shown are from the model calculations of Bennhold and Mart (BM) [45] (dot-dashed), Janssen et al. (JB) [46'] (solid), and Guidal et 
al. (GLV) [47] (dashed). The models are described in the text. The data (and calculations) for ajj for the two back-angle points have been 
scaled by a factor of 3 for clarity. 
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curves shown are from the model calculations of Bennhold and Mart (BM) 
al. (GLV) [47] (dashed). The models are described in the text. 
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K^K and final states, in particular, the Pii(1710) and Pi3(1720). The effect of these 

resonances can be seen in the hadrodynamic model calculations (BM model - dot-dashed 
curves, JB model - solid curves), though clearly their strengths at the measured are not 
correct. 

The double-peaking of au for K production at backward 6]^ angles as seen in Fi g . [T6l 



corroborates a similar structure seen in recent photoproduction results [18|, Il9|, |20|, |21 |. 
Within existing hadrodynamic models, the structure just above the threshold region is typ- 
ically accounted for by the known 5*11(1650), Pii(1710), and Pi3(1720) nucleon resonances. 
However there is no consensus as to the origin of the bump feature at ~1.9 GeV that was 

n 

first seen in the K^A photoproduction data from SAPHIR [ISy. It is tempting to specu- 
late that this is evidence for a previously "missing", negative-parity J = 3/2 resonance at 
1.96 GeV predicted in the quark model of Capstick and Roberts [2j. This explanation was 
put forward in the work of Bennhold and Mart 26|], in which they postulated the existence 
of a Di3 state at 1.9 GeV. However, other groups have shown that the same data can also be 
explained by accounting for w-channel hyperon exchanges |50i] or with an additional P-wave 



resonance 



]. From our data, the W spectra of the interference terms, (Ttt and (Tlt-, show 
no clear structures in the region about 1.9 GeV, whereas an s-channel resonance would likely 
be reflected in the structure of the interference terms, particularly (Jtt- Note that the BM 



and JB models include a -Di3(1895) resonant state whose coup 



18, 



ing strength was determined 



19| . Clearly the differences 



from fits to the SAPHIR K^K total cross section data 
between both models and our data indicate that either the resonance parameters are not 
accurate, that more resonant terms are required, or that the bump at 1.9 GeV in our W 
spectra has a non-resonant origin. We conclude that the W dependence of i^+A production 
provides suggestive evidence for baryon resonance activity within the reaction mechanism, 
but that the data in comparison to present models does not allow any simple statement to 
be made. 



In the K^TP channe 



toproduction result [19 



, ajj is peaked at about 1.85 GeV, which also matches the pho- 



20l |. In addition, (Jtt shows a broad feature in this same region. 



These features are consistent with a predominantly s-channel production mechanism. In 
this region, beyond the specific A^* resonances believed to contribute to K^K production 
(and hence are strong candidates to contribute to K'^TP production), there are a number 



of known A* resonances near 1.9 GeV 



51| that can contribute to the K^TP final state, 
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particularly the A*(1900) and A*(1910). These A* states are forbidden to couple to the 
K^K state due to isospin conservation. Current hadrodynamic models seem to indicate 
that both A^* and A* states (see Table [T]) are necessary to describe the existing photo- and 
electroproduction data. 

The comparison of the hadrodynamic model calculations to the data clearly indicates 
that significant new constraints on the model parameters will be brought about when these 
new data are included in the fits. The models do not reproduce au, cttt, ot aiT at any level, 
especially for the K^Tf data. The Regge model tends to underpredict the strength of au 
across the full angular range, which is suggestive of s-channel contributions to this reaction. 
Again the trends of arr and are reasonably well reproduced with the inclusion of only 
the K and K* Regge trajectories. 

3. Dependence 

The data shown in Figs. [16] and [17] were obtained from our 2.567 GeV data set at 
Q^=0.65 GeV^. Our data set at 4 GeV provides a much larger reach and it is in- 
structive to study the W spectra for increasing values of Q^. These data are shown in 
Figs. [T8] and [T9] for one of our backward-angle points (cos 6'^=-0.25) and a more forward- 
angle point (cos ^^=0.35). The interference structure functions (not shown) do not have a 
strong dependence, while ajj shows a smooth fall-off. Note that at (5^=1.0 GeV^, our 
back-angle data do not show the double-peaked structure that was evident at Q^=0.65 GeV^ 
(see Fig. [T6l) and also in our Q^=l.{] GeV^ data at 2.567 GeV (not shown). This could be 
due to our increased W bin width at 4 GeV (100 MeV compared to 50 MeV at 2.567 GeV), 
or could imply a strong e dependence to the resonance strength. 

None of the models reproduces the K^K data in detail. Both hadrodynamic models are 
very poor matches to these data, while the GLV Regge model tends to underpredict the 
strength in our more forward angle point, although it is in fair agreement with the data in 
our more backward angle point. For the K^TP data none of the models shown reproduce 
even the qualitative aspects of the data. 

The dependence of ajj for K^K and K^TP can be studied within our 4 GeV data 
set as shown in Fig. [20] The data shown are from our points at cos6'j^=-0.25 and 0.90 for 
three different W values across the nucleon resonance region. Also included on these plots 
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FIG. 18: W distributions of au for the K^A final state from our 4 GeV data set (evolved to 
4.056 GeV) for each of our four points in for cos ^^=-0.25 (left) and 0.35 (right). The relative 
systematic uncertainties to ajj are given in Table llVi The curves shown are from the model 
calculations of Bennhold and Mart (BM) [45|] (dot-dashed), Janssen et al. (JB) [46,] (solid), and 



Guidal et al. (GLV) (dashed). 

are our two data points from the 2.567 GeV data set at Q^=0.65 GeV^ and 1.00 GeV^, 
along with the CLAS cr-r data from photoproduction at Q^=0 from Bradford et al. [21]. No 
clear features are apparent here, with the data showing a smooth fall-off with respect to the 
photon point with increasing for both final states except for the forward-angle K^A data 
at W=1.75 GeV. In order to compare more directly with the existing measurements from the 
1970's (taken for 9}^ < 15°, 0.5 < < 4.0 GeV^, and evolved to W=2.15 GeV) compiled by 
Bebek et al. in Ref. jlll, we have fit our 4 GeV au data with the dipole form C/(Q^ + M^)^ 
(where C is an arbitrary constant and the CLAS photoproduction data are not included in 
the fits) and compared the mass terms to those extracted from the fits in Ref. ll|. Fits 
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FIG. 19: W distributions of au for the K^Y? final state from our 4 GeV data set (evolved to 
4.056 GeV) for each of our four points in for cos0^=-O.25 (left) and 0.35 (right). The relative 
systematic uncertainties to ajj are given in Table llVi The curves shown are from the model 

45( 1 (dot-dashed), Janssen et al. (JB) (solid), and 



calculations of Bennhold and Mart (BM) 
Guidal et al. (GLV) (dashed). 

to the older data suggested that the data with = (0.785 ± 0.095) GeV^ fell off 

more rapidly with increasing than the K^A data with = (2.67 ± 0.28) GeV^. The 
results from our fits are contained in Table |V] and shown in Fig. [201 Our extracted mass 
terms, even for backward angles where s-channel and w-channel contributions are expected 
to be more important relative to t-channel kaon exchange, are consistent with the fits of 
Ref. [ul extracted from forward kaon angle data. These results highlight the fact that the 
production mechanisms for K~^A and K~^J]^ are quite different. 

It is interesting to see that the fits for the K'^A data significantly overshoot the photon 
point for our forward angle data. In the absence of other knowledge, one might speculate 
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TABLE V: Mass terms from the fit to our ajj structure functions vs. (not including the 
photoproduction points). A dipole form of C{Q'^ + M^)~^ is employed. 
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FIG. 20: (Color online) distributions of au for the K and K^TP final states from our 4 GeV 
data set (dark filled circles ~ evolved to 4.056 GeV) at cos6l^=-0.25 (left) and 0.90 (right) for 
iy=1.750 GeV (top), 1.850 GeV (middle), and 1.950 GeV (bottom). The solid curves are from a 
dipole mass fit to the 4 GeV data of the form C{Q'^ + M^y^. The Q^=0 points (solid squares) 



come from Bradford et al. 



2l[ and the two data points from our 2.567 GeV data (light crosses) 



are not included in the fits. The dashed lines represent the error bands from the fits. 
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that this is entirely due to a significant contribution to the cross section from cr^. However, 
when the points from the 2.567 GeV data set are included on the plot, we see that they fall 
near the curve fit to our 4 GeV data for all points for K^K and K^TP except our highest W 
point for K^K. This suggests a small contribution from ai. Indeed, the ai/cFT separations 



shown in the next section verify this. Janssen et al. [251] have calculated a fall off in ctt near 
Q^=0, as indicated by our data. This was accomplished by including a dependence to 
the kaon and proton form factors. 

B. ut and Separation 

Our analysis results for ax and ai are presented in Fig. [2T] for the K^A final state and in 
Fig. [21 for the final state. The data are shown here in terms of the ratio R = a^/ (Jt 

as a function of cos^|^ for our different W values. For the K^K final state our analysis 
includes W points from 1.65 to 1.95 GeV, and for the K^TP final state our analysis includes 
W points from 1.75 to 1.95 GeV. Note that the statistical quality of our data did not allow 
us to separate and ax at iy=2.05 GeV. The figures show the ratio extraction using both 
the Rosenbluth and the simultaneous e — $ fit techniques, and the error bars show both 
statistical and total statistical and systematic uncertainties. The discussion of systematic 
uncertainties on these quantities is included in Section IVIII 

The agreement between the Rosenbluth and simultaneous e — $ fits is generally very good 
across our full W and cos6'|^ phase space at = 1.0 GeV^. The ratio of ctl/ctt for both 
the K^A and final states shows ctl to be consistent with zero over our full kinematic 

range except in our highest W point for the K^A reaction, where the value of R varies 
between 0.5 and 1 depending on kaon angle. While several of the extracted values for R are 
negative and might be considered "unphysical" , the majority of these points are consistent 
with zero within the combined statistical and systematic uncertainties. 

Our data at W=1.85 GeV and cos ^^^=0. 90 are consistent with the parallel kinematics 
measurement of Mohring et al. [lj| from Hall C which found R = 0.45lg}g for K^A and 
R = 0.29^|]|!_for K+T,^ at 1^=1.84 GeV, as well as the recent Hall B results of Raue 



and Carman [jj for K+A which found i?=0.005±0.228, 0.239±0.343, and 0.088±0.480 for 
9}^=0° and W=1.72, 1.84, and 1.98 GeV, respectively. Note that the quoted uncertainties 
on R from Mohring et al. [L^ and Raue and Carman given here represent the total 
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statistical and systematic uncertainties. 
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461] (solid), and Guidal et al. (GLV) (dashed). The parallel 



kinematics data points come from Mohring et al. ^14] (open square) and Raue-Carman 
triangles). 



171] (open 



The predictions of the hadrodynamic models for the ratio R are very sensitive to the 
dynamics included in the calculation. All models shown for the i^+A final state predict R 
to be less than unity. Given the size of the error bars on the data, all of the models can 



be said to be roughly consistent with the data. However it is clear that the JB model [46] 
predicts too much longitudinal strength at low W . Also for the highest W point, where R 
begins to increase, the BM model [451] predicts too little longitudinal strength. The GLV 



model 



471 ] is in very good agreement with the data over the full kinematic range shown. 



For the K'^TP final state, the JB model 



461 ] disagrees with the measured ratio i?, which 
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FIG. 22: Results for the ratio R = glIot for the K^Y? reaction for the Rosenbluth technique 
(diamonds) and the simultaneous e — <I> fit (filled circles). The Rosenbluth results have been offset 
in angle for clarity. The data are plotted versus cos 6*1^ for our three W points at (5^=1.0 GeV^. 
The inner error bars are statistical only and the outer error bars are combined statistical and 
systematic. The curves shown are from the model calculations of Bennhold and Mart (BM) 45 1 



(dot-dashed), Janssen et al. (JB) [46\ (solid), and Guidal et al. (GLV) (dashed). The parallel 
kinematics data point comes from Mohring et al. [14] (open square). 



will be shown to result from too little transverse strength. The BM model 45| also suffers 
from too little transverse strength for our two lowest W points, but is consistent with the 
data at W=1.95 GeV. Again, the GLV model [47| agrees well with the data over the full 
kinematic range shown. 

Clearly these data, even with their sizeable statistical and systematic uncertainties, can 
provide for important constraints on the underlying dynamics and production models for 
both the K^A and K^T,^ final states. 

The structure functions ax and ctl are plotted separately in Figs. [23] and [21] for K~^A 
as a function of cos^^^. Here ar has a similar trend in angle for all W points, peaking at 
forward kaon angles and falling off smoothly as the angle increases. For ai, the strength 
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is consistent with zero for our points in W from 1.65 to 1.85 GeV. At W^=1.95 GeV, is 
comparable with ar in its angular dependence and its strength. The data for ctt and ai are 
consistent with the existing parallel kinematics measurement at W=1.85 GeV of Mohring 
et al. [l^ from Hall C at JLab. 

The comparison of the models to o"t for the K~^A final state shows that they underpredict 
the data and the strength of the forward-angle rise for our two lowest W points. For our two 
highest W points, the calculations (with the exception of the JB model 46|] at W=1.95 GeV) 
are in good agreement with the data. For ctl, the calculations from the BM model and 
the GLV mod el |47l| are in reasonable agreement with the data given the error bars. However, 
the JB model j46| predicts too much longitudinal strength for the full kinematic range shown. 

400 
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FIG. 23: Structure function ax vs. cosO'^ for the K^A final state for our different W points at 
Q^=1.0 GeV^ from the e — ^ fit. The inner error bars are statistical only and the outer error 
bars are combined statistical and systematic. The curves shown are from the model calculations 
of Bennhold and Mart (BM) l^sl (dot-dashed), Janssen et al. (JB) 



461 (solid), and Guidal et al. 

n 

(GLV) (dashed). The parallel kinematics data point comes from Mohring et al. [IJ] (open square). 



The structure functions ax and 0"^ are plotted separately in Figs. [23 and [22] for as 
a function of cos^^. cxt is seen to have a broad peaking at more central angles for i^+S^ 
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FIG. 24: Structure function cjl vs. cos9*^ for the K^K final state for our different W points at 
(5^=1.0 GeV^ from the e — $ fit. The inner error bars are statistical only and the outer error 
bars are combined statistical and systematic. The curves shown are from the model calculations 

45( 1 (dot-dashed), Janssen et al. (JB) 



of Bennhold and Mart (BM) 



461 1 (solid), and Guidal et al. 

n 



(GLV) (dashed). The parallel kinematics data point comes from Mohring et al. [IJ] (open square). 

compared to the K final state and ctl is consistent with zero everywhere. The data for ctt 
and ai are consistent with the existing parallel kinematics measurement at PF=1.85 GeV 
of Mohring et al. from Hall C at JLab. 

The models compare poorly with the data for 0"^ for K^TP underpredicting the strength 
of the data and missing the trends in the angular dependence. For ctl, all models predict a 
small strength in agreement with the data. Given the size of the statistical and systematic 
error bars on the data, not much more can be said with respect to the model predictions. 



IX. CONCLUSIONS 

We have measured K^K and K^TP electroproduction over a wide range of kinematics 
in the nucleon resonance region, including first-reported measurements over the full range 
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cos 6^ 

FIG. 25: Structure function vs. cos0|^ for the K^Tj^ final state for our different W points 
at Q^=1.0 GeV^ from the e — $ fit. The inner error bars are statistical only and the outer error 
bars are combined statistical and systematic. The curves shown are from the model calculations 
of Bennhold and Mart (BM) 



45( 1 (dot-dashed), Janssen et al. (JB) 



461 (solid), and Guidal et al. 

n 

(GLV) (dashed). The parallel kinematics data point comes from Mohring et al. [IJ] (open square). 



of kaon angle. We have presented data for the separated structure functions err, ctl, and 
the interference structure functions (Tlt and cttt- We conclude that K^K and K^TP elec- 
troproduction dynamics are markedly different. We find that ajj (= + ecr^) and (Jtt are 
forward-peaked (in kaon angle) for the K^K final state and peaked at more central angles 
for K^Tf with cttt a significant fraction of ctjj for both hyperon final states. For the K^A. 
channel, aiT is a significant fraction of au and negative, while it is consistent with zero for 
the K^TP channel. 

The W dependence of ajj for K^TP shows a broad enhancement at ly ~ 1.85 GeV 
for (Tu and (Ttti presumably due to the various A* resonances in this mass region. The 
W dependence for K^K production is more complicated, evolving from a single-peaked 
structure at forward angles to a double-peaked structure at backward angles at low Q^, with 
the double-peaking not obvious above (5^=1.0 GeV^. The longitudinal structure function 
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FIG. 26: Structure function vs. cos0|^ for the K~^T,^ final state for our different W points 
at Q^=1.0 GeV^ from the e — $ fit. The inner error bars are statistical only and the outer error 
bars are combined statistical and systematic. The curves shown are from the model calculations 
of Bennhold and Mart (BM) 



45( 1 (dot-dashed), Janssen et al. (JB) 



461 (solid), and Guidal et al. 

n 

(GLV) (dashed). The parallel kinematics data point comes from Mohring et al. [IJ] (open square). 

aL is consistent with zero for K^K across our full kinematic space except for our highest 
W point at 1.95 GeV, and is consistent with zero everywhere for our K^TP data. The 
transverse structure function (Jt is forward-peaked for K^K and peaked at more central 
angles for K^TP . 

The dependence of the structure functions is unremarkable; the relatively slow fall-off, 
presumably reflecting the form factors of the various exchanged kaons in the case of t-channel 
processes, and the analogous baryonic form-factors in s and w-channel processes. Of interest 
is our observation that an extrapolation of au to = overshoots the photoproduction 
value for aT- The obvious conjecture that this reveals the presence of a large value of is 
not consistent with our direct measurements of 0"^, at least for lower values of W . 

Detailed calculations are needed in order to investigate whether our data indicate the 
presence of any new N* resonances in K^K production. At the same time, the calculations 
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must be able to fit the strong peak seen for X+E^ production and tfie sfiape of utt and 
a^T in both channels. The question of the presence of any new resonances must wait for 
further work with the existing hadrodynamic models and partial wave analyses applied to 
the full range of our data. Fortunately, the new information we present here, especially 
the interference terms, will impose stringent constraints on the amplitudes used to model 
electroproduction of K^K and K'^Tf final states, making these models more reliable for 
future interpretation and prediction. 
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